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SUMMARY 

A  theoretical  investigation  is  undertaken  to  study  the  buckling  of  a 
thin,  rectangular,  homogeneous,  isotropic  flat  plate  augmented  with 
a  centrally  located  longitudinal  stiffener.   The  influence  on  overall 
plate  efficiency  of  the  use  of  advanced  composite  materials  in  the 
stiffener  is  developed  and  discussed.   Results  of  the  analysis  are 
presented  in  the  form  of  charts  to  demonstrate  the  effects  of  fibrous- 
composite  reinforcement  of  the  stiffener  on  the  critical-stress-to- 
weight  ratio  of  the  plate. 
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INTRODUCTION 

The  trend  to  higher  performance  aircraft  in  all  speed  regimes  has  led 
designers  to  seek  methods  of  decreasing  the  structural  weight  of  air- 
craft while  maintaining  or  increasing  the  stress  levels  on  the  mate- 
rials involved.   One  of  the  most  promising  means  of  accomplishing  this 
objective  has  been  the  introduction  of  advanced  composite  materials  to 
either  augment  or  replace  conventional  structural  materials  (for  exam- 
ple, aluminum).   One  type  of  composite  material  consists  of  the  in- 
corporation of  high-strength,  high-modulus  filaments  into  a  compatible 
matrix.   Several  material  system  classes  currently  of  active  interest 
include  graphite  epoxy,  boron  epoxy,  and  metal  matrices  having  metal 
fibers.   As  pointed  out  in  Reference  1,  design  and  structures  engi- 
neers, using  composites,  have  the  ability  to  build  directional  prop- 
erties into  the  structure  which,  in  turn,  results  in  improved  efficiency, 
Weight  saved,  with  its  direct  effect  on  the  design  index  of  strength/ 
weight,  is  the  most  attractive  feature  of  the  use  of  high-performance 
composites.   It  has  not  been  uncommon  for  structures-materials  engi- 
neers to  predict  weight  reductions  of  as  much  as  forty  percent  in  total 
systems  weight  through  the  potential  use  of  advanced  composites  in  both 
primary  and  secondary  structures. 

With  advanced  filamentary  composites,  increased  design  flexibility  of 
significant  degree  exists  that  is  not  possible  with  conventional  metal 
structures.   However,  as  with  many  new  developments,  work  in  the  com- 
posite field  has  been  primarily  empirical,  while  theory  and  analysis 
have  lagged  somewhat  behind.   Despite  the  potential  of  composites, 
evidence  of  their  use  appears  only  in  the  reasonably  small  structural 
elements  that  can  be  designed  by  test  (for  example,  stabilizers,  wing 
tips,  slats,  and  rotors).   The  trial  and  error  procedure,  especially 
with  materials  of  such  high  performance  possibilities  as  composites, 
can  lead  to  excessive  costs  due  to  both  inherent  high  material  costs 
and  eventual  appearance  of  structural  behavior  neither  anticipated  nor 
properly  estimated.   One  such  case  is  presented  by  Durlofsky  and  Mayers 


in  Reference  2.   Another  such  case  exists  with  the  use  of  advanced 
composites  to  enhance  the  buckling  efficiency  of  stiffened  plates 
utilized  in  the  construction  of  wing,  fuselage,  and  empennage  struc- 
tures of  aircraft.   Peterson  studied  the  effect  of  an  all-composite 
stiffener  in  the  reinforcement  of  a  simply  supported  aluminum  plate 
without  considering  the  possibility  of  an  optimized  composite-metal- 
stiffener  combination.   His  results  indicate  that  some  caution  must 
be  exercised  in  using  high-performance  composites  in  combination  with 
conventional  structural  concepts  and  prediction  methods. 

Earlier  studies  of  conventional  stiffened  plates  under  compression, 
such  as  those  of  Barbre,   Seide  and  Stein,   and  Gomza  and  Stein, 
neglect  the  effects  of  torsional  stiffness  of  the  stiffener  in  es- 
tablishing buckling  criteria.   The  justification  offered  is  usually 
that  the  torsional  stiffness  of  the  stiffeners  is  relatively  small 
when  compared  with  the  bending  stiffness  and  that  vanishingly  small 
torsional  stiffness  of  the  stiffeners  is  consistent  with  the  appear- 
ance of  buckle  nodes  at  stiffener  locations  in  the  chordwise  direc- 
tion.  Such  justification,  of  course,  is  appropriate  for  the  discrete 
stiffener  spacing  of  low  subsonic  speed  flight  vehicles.   For  the 
higher  speed  ranges,  stiffeners  are  more  closely  spaced  and  the 
stiffener  torsional  stiffness  becomes  significant  in  buckling  anal- 
ysis.  Thus,  the  introduction  of  composite  materials  into  all  or  part 
of  a  stiffener  can  possess  ramifications  not  evident  in  the  use  of 
conventional  materials.   For  example,  although  Young's  modulus   (E) 
for  a  composite  may  be  as  much  as  four  times  greater  than  that  of  a 
structural  metal,  the  shear  modulus   (G)   may  be  lower  by  approxi- 
mately the  same  factor.   Therefore,  neglecting  torsional  rigidity 
in  a  composite-metal  system  with  closely  spaced  stiffeners  may  in- 
troduce unsatisfactory  errors  in  buckling  calculations. 

It  is  the  purpose  of  this  investigation  to  (l)  present  the  theory  and 
an  approximate  analysis  for  determining  the  stability  of  a  simply  sup- 
ported, flat,  rectangular  orthotropic  plate  with  a  centrally  located, 


composite-augmented,  longitudinal  stiffener  and  (2)  determine  for  a 
square  plate  the  effects  of  varying  the  percentage  (by  volume)  of 
composite  material  in  the  stiffener  with  stiffener  torsional  rigidity 
included. 


GENERAL  THEORY 

STRAIN-DISPLACEMENT  RELATIONS 

The  strain-displacement  relations  for  small  strain,  large  displacement, 
and  moderate  rotation  behavior  of  an  ortho tropic  plate  (see  Figure  l) 
can  be  obtained  from  Reference  7  as 
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With  a  single  stiffener  buckling  into  no  more  than  2  half-wavelengths 
and  a  homogeneous,  isotropic  base  plate,  however,  interlaminar  and 
transverse  shear  effects  become  of  higher-order  significance  and  they 
are  assumed  to  be  negligible  throughout  the  analysis.   Thus,  the  trans- 
verse shear  angles  j        and  y        are  set  equal  to  zero.   Further,  since 
buckling  is  of  interest,  the  stretching  of  the  middle  surface  of  the 
plate  (and  stiffener)  may  be  neglected;  thus 
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STRESS -STRAIN  RELATIONS 

For  an  orthotropic  plate  in  which  the  stresses  do  not  exceed  the  elas- 
tic limit,  Hooke's  law  is  applicable;  it  is  expressed  in  the  following 


form  relating  normal  force  per  unit  length  of  cross  section  to  strain: 
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A  single  longitudinal  stiffener  possesses  the  same  stress-strain  re- 
lationships as  a  one-dimensional  beam;  that  is 
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POTENTIAL  ENERGY 

The  strain  energy  stored  in  an  elastic  plate  of  constant  thickness  due 
to  bending  and  twisting  may  be  written  as 
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Substitution   of   Equations    (2)    and    (3)    into    Equation    (5)    and    subsequent 
integration  over   the    thickness   of   the   plate  yields 


U 


L        b 

Iff 


0        0 


D  _             D  [i      +  D  n 
x  2               x^y          y  x 
w  +  * * — 

i  }  xx  . 

1-UU  1     -     U    (.1 

x^y  xh  y 


w        w 
,xx   ,yy 


(Continued) 


1  -  n  |i 

x  y 


2  2 

w    +  2D  w 
yy     xy  ,xy 


dxdy 


(6) 


With  the  use  of  the  Dirac  delta  function,  the  potential  energy  of  a 
single  centrally  located  symmetric  longitudinal  stiffener  (see  Fig- 
ure l)  may  be  described  over  the  same  interval  as 
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Substitution  of  Equations  (2)  and  (k)    in  Equation  (7)  and  subsequent 
integration  over  the  stiffener  height  yields 
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If  an  axial,  inplane  load,  positive  in  compression,  is  applied  to  the 
edges  of  the  plate,  the  potential  of  this  applied  load  (including  the 
end  load  on  the  stiffener)  is 
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However,  since  in  the  absence  of  middle  surface  stretching  (e   =  e 

>  \         X  y 

7    =0  when   z  =  0),  it  can  be  observed  from  Equation  (l)  that 
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Substitution  of  Equation  (10)  in  Equation  (9)  yields 
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Finally,  the  total  potential  is  the  sum  of  Equations  (6),  (8) ,  and  (ll) 
and  becomes 
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VARIATIONAL  PRINCIPLE  AND  GOVERNING  EQUATIONS 


In  accordance  with  the  potential  energy  principle,  the  total  potential, 
U  +  V   ,  must  have  a  stationary  value  with  respect  to  admissible  var- 
iations in  the  lateral  deflection,   w   ,  in  order  for  equilibrium  to 
be  achieved.   This  variation  leads  to  the  equilibrium  or  Euler  equation 


of  the  variation,  plus  the  applicable  boundary  conditions.   The  var- 
iation of  Equation  (12)  leads  to  the  following  equilibrium  equation: 
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The  associated  boundary  conditions  for  the  case  of  single  support  are 
obtained  as 
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The  solution  of  Equation  (13);  consistent  with  the  boundary  conditions, 
Equation  (lU)  leads  to  the  deflection  function,   w   .   This,  however, 
is  a  difficult  procedure  in  general  and  is  not  the  approach  followed 
in  this  investigation.   Rather,  the  direct  method  of  the  calculus  of 
variations  (Rayleigh-Ritz  method)  is  used  to  determine  the  buckling 
criteria  for  w  expressed  in  terms  of  a  series  of  arbitrary  parameters 


METHOD  OF  SOLUTION 

With  the  total  potential  energy  given  by  Equation  (12)  in  terms  of  the 
plate  deflection  w   ,  the  direct  method  of  the  calculus  of  variations 
(Rayleigh-Ritz  method)  is  used  to  determine  the  minimum  buckling  load 
for  varying  stiffener  size  and  composite  content. 

GENERAL  CASE 

A  deflection  function  satisfying  the  boundary  conditions  for  a  flat, 
simply  supported,  stiffened  rectangular  plate  under  axial  compression 
is 


w 


Z.  .      mrrx      .      nrry 

A        sin  — —  sin  — r-*- 
mn  L  b 


(15) 


n 


Substitution  of  Equation  (15)  into  Equation  (12)  (with  N  =  N   =  0  ) 
and  subsequent  integration  over  the  surface  of  the  plate  yields  (see 
Appendix  i) 
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The  vanishing  of  the  first  variation  of  the  total  potential  energy  with 
respect  to  A    ,  the  nondimensional  arbitrary  displacement  parameter, 
yields  the  equilibrium  condition  and  a  criterion  for  obtaining  the 
buckling  load.   This  procedure,  as  shown  in  Appendix  I,  yields 
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HOMOGENEOUS,  ISOTROPIC  SQUARE  PLATE 

A  homogeneous,  isotropic  square  plate  with  an  equivalent  stiffener  (to 
account  for  varying  amounts  of  composite  material)  has  been  selected 
as  the  specific  case  for  examination  in  the  present  investigation.   A 
square  plate  with  one  centrally  located  longitudinal  stiffener  under 
these  loading  conditions  would  buckle  in  one  of  two  configurations: 
(l)  one  half-wave  in  each  direction   (m  =  n  =  l)   or,  (2)  two  half- 
waves  in  each  direction   (m  =  n  =  2)   .   Any  other  combinations  of 
m  and   n  would  not  lead  to  the  minimum  value  of  the  buckling  co- 
efficient k    .   As  noted  from  Equation  (IT)  the  bending  stiffness 
term   (El)/Db   of  the  stiffener  is  effective  only  during  the   1,1 
mode,  and  conversely,  the  torsional  stiffness  term   (GJ)/Db   is  ef- 
fective only  during  the  2,2  mode. 

For  an  approximate  solution,  the  first  term  of  the  deflection  function, 
Equation  (15);  could  be  used.   However,  this  would  allow  only  for  a  pure 
sinusoidal  half-wave  deflection  pattern  of  the  plate  in  the  chordwise 
direction.   As  stiffener  size  or  stiffness  is  increased,  rotation  at  the 
inflection  point  of  the  2,2  mode  is  somewhat  restrained,  indicating  the 
need  for  at  least  a  two-term  deflection  function;  that  is,   n  =  1,3   for 
the   1,1  mode  and   n  =  2,k      for  the  2,2  mode  of  Equation  (17).   The 
result  of  carrying  out  this  procedure  (see  Appendix  I)  is 
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for  the   1,1  mode  and 
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for  the  2,2  mode,  where   D   is  the  bending  rigidity  for  a  homo- 
geneous, isotropic  plate.   The  quantities   (El)   and   (GJ)   are  the 
bending  and  torsional  stiffnesses,  respectively,  of  the  equivalent 
longitudinal  stiffener.   The  procedure  for  determining  these  values 
for  varying  amounts  of  composite  material  in  the  overall  stiffener 
is  given  in  Appendix  II. 

The  minimum  values  of   k   may  be  calculated  from  Equations  (l8)  and 

(19)  in  each  of  the  two  modes  of  the  square,  stiffened  plate  for 

varying  values  of   (El)/Db  ,  (GJ)/Db  ,  and   (A  )/bt   .   The  critical 

stress,   cr    ,  is  then  calculated  from  the  relationship 
'        cr   ;  r 


cr 


2 
k  7T  D 
x 


b2t 


(20) 


The  average  density  of  the  plate-stif fener  system  is  a  nondimensional 
measure  of  the  relative  weight  of  the  combination  with  the  introduc- 
tion of  composite  materials  into  the  system.   This  factor  is  calculated 
by  averaging  the  densities  of  the  system  as  follows 

p     =   p     K   comp)  +  p,     (<fo   base)  (21) 

ave      comp  u     r/    base  w      '  v  ' 

where   ($  comp)   is  the  percent  by  volume  of  composite  in  the  entire 
plate  stiffener  system,  and   ($>  base)   is  the  percent  by  volume  of  the 
base  metal  used  (aluminum,  steel,  etc.). 

The  buckling  efficiency  of  the  stiffened  plate  combination   (cr  /p    ) 

to  J  r  \  cr/  ^ave  ' 

is  plotted  versus   (A  )/bt   for  varying  amounts  (by  volume)  of  composite 
material  in  the  stiffener  in  Figures  2  through  k.      Aluminum  is  the  base 
material,  with  either  boron  epoxy  or  graphite  epoxy  used  as  the  advanced 
composite  material. 

Peterson  determined  that  the  optimum  longitudinal  stiffener  for  a  flat 
plate  has  the  relative  proportions  (l)  stiffener  width  equal  to  twice 
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plate  thickness,  and  (2)  stiffener  height  equal  to  eight  times  the 
plate  thickness.  These  ratios  are  utilized  as  starting  values  for 
the  present  investigation. 
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RESULTS  AND  DISCUSSION 

An  investigation  of  the  axial  compression  buckling  of  a  simply- 
supported  flat  plate  with  a  central  longitudinal  stiffener  has  been 
undertaken  to  estimate  the  effects  of  composite  reinforcement  of  the 
stiffener  with  stiffener  torsional  stiffness  effects  included.   The 
minimum  or  critical  values  of   k    ,  using  Equations  (l8)  and  (19);  are 
calculated  for  the  two  modes  of  a  stiffened,  square,  homogeneous,  iso- 
tropic plate  for  varying  values  of   (El)/Db   and  plotted  in  Figure  5« 
These  values  are  in  excellent  agreement  with  those  of  Seide  and  Stein 

for  the  m  -   1  ,  n  =  1   case,  but  the   (  xj    values  of  this  investi- 

7  '  x   'cr 

gation  are  slightly  higher  than  those  of  Reference  5  for  the  m  =  2  , 

n  =  2   case.   This  discrepancy  is  due  to  the  inclusion  of  the  torsional 

stiffness  of  the  stiffener  in  the  present  analysis.   Without  torsional 

stiffness  effects,   I  x)    for  the  m  =  2  ,  n  =  2   case  remains  con- 
>      \       /cr  > 

stant  at  a  value  of   16  despite  the  size  of  the  stiffener.   Timoshenko 

o 
and  Gere,   in  the  case  of  a  plate  simply  supported  on  three  edges  with 

one  of  the  unloaded  edges  clamped  under  compressive  loading  show  the 

maximum  value  of  \    xj     to  be   22.6  .   This  would  be  analogous  to  the 

present  problem  when  the  torsional  stiffness  is  infinitely  large.   The 

two-term  deflection  function  used  in  the  present  investigation  will 

result  in  values  of  \    xj  higher  than  that  arrived  at  by  Reference  8, 

but  only  when  the  stiffener  size  is  excessively  large  and  out  of  the 

area  of  practical  interest.   Equations  (l8)  and  (19)  give  values  of 

\    xj    which  are  sufficiently  accurate,  as  the  primary  aim  of  the 

present  report  is  to  investigate  the  trends  of  plate  efficiency  rather 

than  to  calculate  actual  design  data  concerning  buckling  loads. 

The  primary  quantity  for  demonstrating  plate  buckling  efficiency  in  this 

investigation  is   cr   /p      ,  where   cr    is  the  critical  stress  of  the 
b  cr'  Kave   '         cr 

plate  determined  from  Equation  (20),  and   p      is  the  average  density 
determined  from  Equation  (21).   Figure  2  is  a  graph  of   cr  /p     for  a 
square  aluminum  plate  with  a  centrally  located  longitudinal  stiffener 
augmented  by  varying  amounts  of  graphite-epoxy  composite.   The  buckling 
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efficiency  has  been  calculated  for  varying  stiffener  size.   The  sharp 
knee  in  each  of  the  curves  results  from  the  transition  of  the  plate 
buckling  pattern  from  a  m  =  1  ,  n  =  1   configuration  to  a  m  =  2  , 
n  =  2  mode.   To  the  left  of  this  knee,  or  in  the   1,1   mode,  sig- 
nificant increases  in  buckling  efficiency  result  from  the  replacement 
of  aluminum  with  graphite  epoxy  in  the  stiffener.   This  is  due  to  the 
increased  flexural  stiffness  of  the  stiffener,  since  the  graphite  epoxy 
has  a  greater  value  of  Young's  modulus  than  does  aluminum.   However,  to 
the  right  side  of  the  knee,  in  the  2,2  buckling  configuration,  the 
plate  buckling  efficiency  is  actually  reduced  with  the  addition  of 
the  graphite  epoxy.   Once  the  2,2   condition  is  reached,  the  stif- 
fener lies  at  a  node  of  the  plate,  and  the  flexural  stiffness  of  the 
stiffener  has  no  effect  on  the  magnitude  of  the  critical  stress.   Only 
the  torsional  stiffness  of  the  stiffener  is  a  factor  in  this  range  and, 
since  the  shear  modulus  of  aluminum  is  approximately  three  times  that 
of  graphite  epoxy,  the  overall  buckling  efficiency  decreases  with  the 
addition  of  the  composite  material.   As  an  example,  a  plate  with  an 

all-aluminum  stiffener  and  a  value  of   (A  )/bt   of   0.205   attains  a 

s        i 

buckling  efficiency  parameter  value  of  ^1.0  x  10   in.  (Point  A).   If 
the  aluminum  stiffener  is  then  replaced  by  an  all  composite  stiffener, 
this  parameter  value  is  decreased  to  ^>Q.Q   x  10   in.  (Point  B). 

It  may  also  be  noted  from  Figure  2  that  in  the   1,1   configuration,  as 

the  percentage  of  composite  in  the  stiffener  increases,  the  percent 

increase  in  the  buckling  efficiency  decreases.   For  a  plate  with  a 

value  of   (A  )/bt   of   0.08  ,   <x  /p     equals   15. 7  x  10   in.   for 
v  sJ '  '        cr  ave 

an  all  aluminum  stiffener.   The  replacement  of  10  percent  of  the  stif- 
fener cross-section  with  graphite  epoxy  raises  the  value  of  cr  /p 

i  to   r  cr  ave 

to  19«5  X  10  in.,  a  percentage  increase  of  almost  25  percent.   Raising 


the  percentage  of  composite  to  50  percent  gives  a  buckling  efficiency 
parameter  value  of  28.3  X  10   in.  ,  or  an  increase  of  80  percent 
over  the  all  aluminum  stiffener.   An  all  composite  stiffener  has  a 

/  x-    h 

cr   /p     value  of  31.6  x  10   ,  or  a  little  greater  than  a  100  per- 
cr  ave  F 

cent  increase  over  the  metal  stiffener.   In  another  case,  where  the 
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size  of  the  stiffener  is  increased  to  give  a  value  of   (A  )/bt  of 
0. 1^  ,  a  10-percent-composite  stiffener  has  approximately  the  same 

value  of  cr  /p     as  the  all-composite  stiffener,  although  the 

cr'  ave  r  '  b 

latter  structure  will  be  lighter  in  weight  and  much  higher  in  cost. 

As  mentioned  previously,  past  studies  having  to  do  with  the  buckling 
properties  of  all-metal  stiffened  plates  have  generally  neglected  the 
torsional  rigidity  of  the  stiffener.   Figure  3  is  a  graph  of  the  same 
plate-stiff ener  combination  as  that  of  Figure  2;  however,  the  shear 
modulus  of  graphite  epoxy  is  assumed  to  be  identical  to  that  of  alu- 
minum  (G  =  ^.0  x  10  lb/in.  )  .   In  this  case,  the  replacement  of  a 
percentage  of  the  aluminum  in  the  stiffener  with  graphite  epoxy  always 

results  in  an  increase  in  cr  /p     without  regard  to  the  percentage 

cr  ave  b  v  b 

of  composite,  size  of  the  stiffener,  or  the  buckling  pattern  of  the 
plate.   Thus,  it  is  evident  that  the  torsional  stiffness  of  a  rein- 
forced metal  stiffener  can  play  a  significant  role  in  the  stability 
problem  and  should  not  be  overlooked  when  high-extensional  modulus 
composites  are  introduced  to  improve  structural  efficiency. 

Another  popular  composite  material  in  developmental  use  is  boron  epoxy. 

This  composite,  although  having  a  higher  value  of  Young's  modulus,  is 

appreciably  more  dense  and  significantly  more  expensive  than  graphite 

epoxy.   Figure  h   is  a  plot  of  the  buckling  efficiency  of  a  boron-epoxy- 

augmented  stiffened  plate.   Comparison  of  Figure  h   with  Figure  2  reveals 

that  although  the  use  of  boron  epoxy  results  in  higher  values  of  cr  /p 

r  °  cr'  ave 

than  graphite  epoxy  in  the   1,1  buckling  configuration,  it  results  in 

greater  loss  in  efficiency  when  in  the  2,2  mode.   Higher  values  of 

cr  /p     in  the   1,1  mode  are  the  result  of  the  increased  Young's 
cr'  ave  '  fe 

modulus  for  boron  epoxy,  whereas  the  reduced  efficiency  in  the  2,2 
configuration  is  attributed  to  the  increased  weight  of  boron  epoxy  and 
the  fact  that  the  shear  moduli  of  the  two  materials  are  essentially  the 
same. 
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CONCLUDING  REMARKS 

This  investigation  develops  a  buckling  criterion  for  axially  compressed, 
simply  supported,  flat,  rectangular,  homogeneous,  isotropic  plates  with 
composite-reinforced,  centrally  located,  longitudinal  stiffeners.   This 
criterion  is  obtained  from  a  more  general  development  based  on  ortho- 
tropic  plate  theory.   Comparison  of  the  results  of  the  present  study 
with  past  investigations  shows  the  derived  buckling  coefficient  to  be 
generally  more  realistic  than  those  which  neglect  the  torsional  stiff- 
ness of  the  stiffener  when  the  plate  cross-section  aspect  ratio  between 
stiffeners  is  at  least  sixty  or  less. 

It  is  also  shown  that  the  addition  of  composites  to  improve  the  effi- 
ciency of  a  stiffened  flat  plate  requires  more  judicious  consideration 
than  simply  replacing  conventional  material  with  high-performance  fi- 
brous composites.   Effects  which  may  have  proved  to  be  negligible  in 
the  analysis  of  conventional  structures  can  become  significant  in  the 
conventional-composite  combination  due  to  the  relatively  large  spread 
between  E   and   G   in  advanced  composites  as  compared  with  that  in 
conventional  metals.   The  replacement  of  a  conventional  stiffener  with 
one  consisting  of  composites  will  usually  lead  to  lighter  structural 
weight,  but  the  undesirable  side  effects  of  lower  buckling  efficiency  and 
increased  cost  may  also  be  introduced.   Careful  study  must  be  made  as  to 
the  amount  of  composite  to  be  used,  the  percentage  of  increase  in  buckling 
efficiency  actually  achieved,  and  the  additional  costs  involved.   It  has 
been  shown  that  in  certain  cases,  a  rather  minimal  use  of  composites  will 
result  in  maximum  gain  in  weight  reduction  and  increase  in  critical  stress 
with  relatively  small  increases  in  cost. 

When  working  with  advanced  composites,  the  designer  of  high-performance 
vehicles  can  no  longer  neglect  the  torsional  rigidity  of  a  stiffener  as 
was  normally  done  with  low  subsonic-range  all-metal  systems.   Such  ne- 
glect can  introduce  errors  which  may  tend  to  show  that  any  use  of  com- 
posites will  achieve  additional  gains  in  structural  efficiency,  when  in 
fact  this  is  not  so.   The  lower  shear  modulus  of  composites  requires  that 

IT 


the  effects  of  torsional  rigidity  should  always  be  taken  into  account,  as 
unconservative  errors  in  critical  stress  may  result  in  either  premature 
buckling  or  reduced  buckling  efficiency. 

The  designer  must  also  be  cognizant  of  all  the  advantages  and  disadvan- 
tages of  the  various  composite  materials  now  becoming  available.   As 
shown  herein  for  graphite  and  boron  epoxy,  each  of  these  materials  has 
certain  advantages,  and  the  most  suitable  material  must  be  chosen  for 
the  specific  application  at  hand.   Composite  materials  present  the  pos- 
sibility of  increased  strength  and  stiffness  along  with  decreases  in 
structural  weight,  but  careful  study  must  be  made  with  regard  to  proper 
use  and  material  selectivity  to  achieve  optimized  structural  efficiency 
and  fabricability  at  minimum  cost. 
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APPENDIX  I 
DEVELOPMENT  OF  BUCKLING  CRITERION 


The  total  potential  energy,  Equation  (12)  (with   N  =  N   =  0  ),  upon 
substitution  of  the  deflection  function,  Equation  (15);  yields 
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Integration  of  this  expression  over  the  area  of  the  plate  yields 
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Now,  the  first  variation  of  Equation  (23)  with  respect  to   A    must 

7  mn 

reduce  to  zero  since  any  first  order  change  in  the  total  potential 
must  vanish  for  every  admissible  variation  in  the  lateral  displacement 
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function.      Hence, 
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Equation    (2U)    reduces    to   Equation    (IT)   when   a   factor   of     m  tt  bD   /klP 
is    removed.      For    the   special    case  of   an   isotropic   square   plate 
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x  y  v 


D  =     D(l   -   \l) 

xy 


(25) 


u        =     u        =     u 
Kx  "y 
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Substitution   of   Equation    (25)    into   Equation    (17)    yields 


mn 


[(l  ♦  n2^ 


k  a 
x 


2(EI)  m    _ 

■  sin  —    2^ 


PIT 


Db 


A        sin 
mp 


n 


2(GJ)     /  b    \  2  nrr     _  /prr  \  pir 

+  [  —    )    n  a      cos  —       /        A       (  —   J    cos  — 

Db         Vim/  2      ^      mp\b    /  2 


2(As) 

k  a     sin 

bt  X 


tL 


pn 

A        sin  —     =      0 
mp  0 


(26) 


A  second   approximation   solution  of  Equation    (26)    for   the  one  half-wave 
pattern   in  each  direction  yields    two   equations    in   two   unknowns:      namely 


n  =  1      :      A  .    (k  -  k  )  + 

ml    ^  x 


2(EI) 


Db 


(A   ,    -  A  _) 

v   ml  m3 


2  (A  ) 

—  k     (A  i    -   A  J 

bt  x       ml  m3y 


=      0 


n   =  3      :      A  2  (100-  k    )  - 
m3  x 


2(EI) 


Db 


(A   ,    -  A  ,  j 
ml  wq 


(27) 


2  (A  J 


bt 


k      (A   .    -    A  ,)     =      0 
x    v   ml  no 
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Now,  for  the  A   's   to  possess  other  than  trivial  values,  the  deter- 
;  mn       r  ; 

minant  formed  by  the  coefficients  of  the   A   's   must  vanish.   Hence, 

mn  ' 


2(EI)   2(A) 
k   -   k  +  -  2- 


Db 


bt 


2(EI)   2(A) 
+  —  k 


Db 


bt 


x 


2(EI)    2(A) 
+  —  k 


Db 


bt 


2(EI)   2(A) 

100  -  k  +  -  k 

x  ,  .     x 


Db 


bt 


=  0 


(28) 


The  solution  of  this  determinant  yields  Equation  (l8).   The  second- 
approximation  solution  of  Equation  (26)  for  the  two  half-wave  pattern 
in  each  direction  yields 


n  =  2 


Am2^ 


2(GJ) 


Db 


(A  -  -  2A  ,,' 
v  m2     m4' 


=   0 


n  =  k 


\k{2k- 


8(GJ)  /a 


m2 


=   0 


Db 


m'- 


(29) 


It  follows  from  the  same  argument  as  before  that  for  a  solution; 


2(GJ)    k 

i    _____    x 
4  +  -  

Db      k 


k(G3) 


k(G3) 


Db 


8(GJ)    k 


25  + 


Db 


Db 


_x 
1+ 


=   0 


(30) 


The  solution  of  Equation  (30 )  yields  the  result  given  in  Equation  (19) 
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APPENDIX  II 

DETERMINATION  OF  EQUIVALENT  VALUES  OF   (El)   AND   (GJ) 
IN  COMPOSITE- REINFORCED  STIFFENER 

o 
Timoshenko  and  Young  describe  a  common  method  of  transforming  a  com- 
posite beam  into  an  equivalent  beam  of  one  material  for  the  determi- 
nation of  bending  stiffness.   This  is  done  by  multiplying  the  width 
of  the  composite  section  of  the  stiffener  by  the  ratio  of  the  Young's 
moduli;  that  is 


E 


w,    ^^  wp  (31) 

E 
base 

where  w   is  the  original  width  of  the  stiffener,  and  wp   is  the 

modified  width  of  the  composite  section  of  the  equivalent  stiffener. 

As  an  example,  suppose  E     =  3E,      .   Then,  as  shown  in  Figure  6, 
r      '         rr      comp     base         '  ' 

the  equivalent  stiffener  would  consist  of  one  material,  the  basic 
metal,  but  it  would  possess  a  different  shape.   The  moment  of  inertia 
of  the  equivalent  stiffener  is  then  found  and  multiplied  by  E,      , 
to  give  the  equivalent  value  of   (El)   for  the  stiffener. 

Timoshenko  and  Goodier   give  a  formula  for  the  twisting  moment  or 
torque  of  a  rectangular  beam  in  torsion  as 

1 


T      -  G0(w  )5  h 
3 


192  w    —^  1  krrh 

1  -  — —  —  y  —   tanh  

T?        h  ,~.  \P  2wn 

k=odd  1J 


(32) 


However,  from  Reference  9 ,     the  torsional  stiffness  of  a  beam  in  free 
torsion  is  known  to  be 

T 

J   =  ~  (33) 

G6 
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Substitution  of  Equation  (32)  into  Equation  (33)  yields 


J   = 


(wx)3  h 


192  w 


7T 


00 
k=odd 


k7Th 


tanh 


2w 


1  J 


(W 


This  value  is  calculated  for  each  section  of  the  stiffener  and  multi- 
plied by  its  respective  shear  modulus.   Addition  of  these  components 
gives  the  equivalent   (GJ)   of  the  stiffener. 
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STIFFENER   AND 
PLATE   X- SECTION 


in 


COMPOSITE 


ALUMINUM- 


Figure  1.   Geometry  of  Homogeneous,  Isotropic  Plate 
With  Longitudinal  Stiffener. 
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Figure  2.   Buckling  Efficiency  Curve  for  Aluminum  Plate  With 
Longitudinal  Stiffener  Augmented  With  Varying  Per- 
centage (by  Volume)  of  Graphite  Epoxy. 
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Figure  3«   Buckling  Efficiency  Curve  for  Aluminum  Plate  With 
Longitudinal  Stiffener  Augmented  With  Varying  Per- 
centage (by  Volume)  of  Graphite  Epoxy  (G     =  G,    ) 
°   x  J  comp    base 
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Figure  k.      Buckling  Efficiency  Curve  for  Aluminum  Plate  With 
Longitudinal  Stiffener  Augmented  With  Varying  Per- 
centage (by  Volume)  of  Boron  Epoxy. 
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